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ABSTRACT 

We present 21 -cm H I observations from the Southern Galactic Plane Survey of the field around the anoma- 
lous X-ray pulsar IE 1048.1-5937, a source whose X-ray properties imply that it is a highly magnetized 
neutron star (a "magnetar"). These data reveal an expanding hydrogen shell, GSH 288.3-0.5-28, centered 
on IE 1048.1-5937, with a diameter of 35 x 23 pc (for a distance of 2.7 kpc) and an expansion velocity of 
« 7.5 km s"' . We interpret GSH 288.3-0.5-28 as a wind bubble blown by a 30^0 M© stai-, but no such central 
star can be readily identified. We suggest that GSH 288.3-0.5-28 is the wind bubble blown by the massive pro- 
genitor of IE 1048.1-5937, and consequently propose that magnetars originate from more massive progenitors 
than do radio pulsars. This may be evidence that the initial spin period of a neutron star is correlated with the 
mass of its progenitor, and implies that the magnetar birth rate is only a small fraction of that for radio pulsars. 
Subject headings: ISM: bubbles — ISM: individual (GSH 288.3-0.5-28) — pulsars: individual (IE 1048.1- 
5937) — radio lines: ISM — stars: winds, outflows — stars: neutron 



\ 1. INTRODUCTION 

The last decade has revealed remarkable diversity in the 
young neutron star population: radio pulsars, soft 7-ray re- 
peaters (SGRS), anomalous X-ray pulsars (AXPs) and cen- 
tral compact objects (CCOs) are all now known to be po- 
\ tenti al compact remnants of core-collapse supernovae (e.g., 
■ iKaspi & Helfand 2002). While the CCOs remain enigmatic, 
the AXPs and SGRs are now both believed to be popula- 
tions of "magnetars", neutron stars whose persistent X-ray 
emission and occasional X- and 7-ray bursts are powered by 
the energy associated with extreme surface magnetic fields, 
> lO'*" G (see (Woods & Thompson 2005. for a review). 

The low Galactic latitudes and the associations of some 
of these sources with supernova remnants (SN Rs) make it 
clear that magnetars are young neutron stars (iGaensler et al.l 
' 1^01). However, what is not clear is why some neutron stars 
are "normal" radio pulsars, while others are X- and 7-ray 
emitting magnetars. The physical distinction cannot simply 
be the strength of the dipole magnetic field, since there is 
now a small population of radio pulsars whose dipole field 
strengths (as inferred from their spin down) are at magnetar- 
like levels, but which have completely different X-ray charac- 
teristics (Pi vovaroff et al. 2000; Kaspi & McLaughlin 2005.). 
Amongst the possibilities suggested to explain the differences 
between these two populations are that the magnetic fields of 
radio pulsars and magnetars differ in their geometry or o ri- 
entation (e.g., 'Zhang & Harding 2000; Kulkarni et al. 2003), 
that magnetars originate from more massive progenitors than 
do radio pulsars (Eikenberrv et al. 2004 ), that high-field radio 
pulsars are quiescent magnetars ( Kaspi & McLaughUnl2005h . 
or that radio pulsars and magnetars for m an evolutionary se- 
quence llLvnei2004HHevl & Hernaiiistl20?)5l) . 
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The environments of magnetars provide unique constraints 
on their origin and evolution. Here we present H I obser- 
vations of the interstellar medium (ISM) surrounding the 
AXP IE 1048.1-5937 (lGavriir& Kaspi 2004, & references 
therein), which show that this source is possibly associated 
with a large expanding H I shell. This association strengthens 
the possibility that magnetars are formed from more massive 
progenitor stars than are radio pulsars. 

2. OBSERVATIONS 

The Hi data shown here were taken from the Souther n 
Galactic Plane Survey (SGPS: lMcClure-Griffiths et alJ2005h . 
which has mapped H I and 1 .4-GHz continuum emission in 
the range 253° < / < 358°, \b\ < l.°5, at a resolution of 
« 2'. The survey combines interferometric data from the Aus- 
tralia Telescope Compact Array with single-dish data from the 
Parkes 64-m radio telescope, giving sensitivity to a wide range 
of spatial scales. The final spectral line data-set has a velocity 
resolution of ~ 0.8 km s"' and a sensitivity of 2 K. 

The SGPS reveals a striking cavity in H I, almost centered 
on the position of IE 1048.1-5937 (namely I = 288.°3,fe = 
-0.°5), as shown in Figure^ The shell has a central velocity 
relative to the Local Standard of Rest (LSR) V = -28 km s"' ; 
we thus designate this object GSH 288.3-0.5-28. At this 
central velocity, the major and minor axes of GSH 288.3- 
0.5-28 are 45' and 29', respectively. The images of / ver- 
sus V and b versus V, also shown in Figure [0 demonstrate 
that GSH 288.3-0.5-28 is expanding, at a velocity V^xp ~ 
7.5 km s"'. GSH 288.3-0.5-28 is seen in multiple velocity 
planes, shows a high contrast between its walls and interior, 
and changes in angular extent as a function of velocity. It thus 
meets all the standard criteria for H I shells in the ISM (see 
pcClure-Griffiths et al. 2002). 

Standard rotation curves (e.g. jBrand & ^Blitz 1993) forbid 
systemic motions having the central velocity observed for 
GSH 288.3-0.5-28, indicating significant deviations from cir- 
cular rotation in this region. However, the shell's central ve- 
locity matches the observed terminal velocity in this direc- 
tion, which allows us to estimate a distance of 2.7 kpc from 
geometry alone. We assume that random cloud motions on 
the order of 6 km s"' contribute to an error in the distance es- 
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timate on the order of 40%, or 1 kpc. Incorporating both this 
distance uncertainty and the non-circular shape of the shell, 
we adopt a radius /? = 14 ± 7 pc. By integrating along a num- 
ber of lines of sight through the shell walls and then subtract- 
ing a mean background, we find an average hydrogen column 
density through the shell walls of (2.5 ±0.6) x 10^" cm"^. 
If this material uniformly filled the region into which the 
shell has expanded, the implied ambient number density is 
Mo~ 17±9cm-3. 

3. INTERPRETATION 

3.1. The Nature ofGSH 288.3-0.5-28 

Shells like GSH 288.3-0.5-28 are common, and result 
from the effect of massive stars on the ISM. Specifically, 
GSH 288.3-0.5-28 could be an H II region, a SNR, or a stellar 
wind bubble. Below we briefly consider these possibilities 

The photoionizing flux needed to maintain an H II region 
is A^* = ATrR^nlaB, where a « 3 x I 0~'^ cm^ s~' is th e re- 
combination coefficient for hydrogen JOsterbrockl 1 9890 . For 
the observed values of R and no, we require A^* ~ (3 ± 2) x 
10^^ s"'. This is consistent with the ionizing fluxes of the 
earliest-type stars, but such a source would also have a radio 
flux of > 100 Jy at decimeter wavelengths. A radio contin- 
uum image of this region sho ws no such em ission associated 
with GSH 288.3-0.5-28 (Gre en etalJl999l) . 

For the observed values of no, R and Vexp, a SNR in the 
Sedov or radiative phases requires an initial explosion energy 
of 5 X 10"*^ or ~ 10^*' ergs, respectively. These estimates 
are both well below the typical supernova explosion energy 
of 10^' ergs, making this interpretation unlikely. There is also 
no evidence in deep archival observations for any radio or X- 
ray SN R associated with IE 1048.1-5937 (see Gaensler et al. 

HooT). 

The remaining possibility is that GSH 288.3-0.5-28 is 
a wind-drive n bubble. Using the wind bubble solution of 
IWeaver et alJ (IT977l) . during the energy-conserving phase of 
expansion we expect that R = Q.16iL„/ pQ^l^t^l^ , where t 
is the age of the bubble, L^. is the wind luminosity of the 
central star, and pQ is the density of the ambient medium. 
Since Vexp = dR/dt, we find L„, w (4±2) x 10^^ ergs s"' and 
t = 1.1 ±0.5 Myr. The luminosity of the wind is a strong 
function of mass, allowing a reasonable typing of the associ- 
ated central object (assuming only one star powers this bub- 
ble). In this case, an appropriate central source is an 06V star, 
with a mass-loss rate M w 2 x 10"^ Mfn yr~' dde Jage r et al.l 
11588'), a wind velocity « 2500 km s"' ('Prinia et al."1990^, and 
a zero-age main sequence ( ZAMS) m ass of 30^0 Mq (e.g., 
iMassev et aU2002l:iOstrov & Lapassel2603i) . 

3.2. GSH 288.3-0.5-28 and Carina OBI 

We have shown that GSH 288.3-0.5-28 is most likely pow- 
ered by the wind of a massive star at a distance of 2.7 ± 1 kpc. 
But can we identify the central star responsible for this bub- 
ble? There are many massive stars in this direction and at this 
distance, most of which are part of th e association Carina OB 1 
at a distance of ^ 2.5 kpc (e.g., Humphreys 1978). However, 
we now argue that GSH 288.3-0.5-28 is unassociated with, 
and probably more distant than, the stars in Car OBI. 

First, the collective winds of the massive stars in Car OB 1 
correspond to a total luminosity of 8 x lO-'^ ergs s"' ( Abbott 
[l282). Combined with the powerful winds from stars in the 
nearby clusters Tr 14 and Tr 16, this should power a much 
larger expanding superbubble. Indeed such a larger shell. 



100-200 pc across, h as been iden t ified in ionized, neutr al and 
molecular gas (Cow ie et al.ll98UlRizzo & Arnall998h . Any 
individual star in this vicinity will not sweep up a neutral shell 
of its own but will contribute to the overall ionized supershell. 
Since the properties of GSH 288.3-0.5-28 are consistent with 
a single star sweeping up neutral gas, it seems unlikely that 
this shell is in the vicinity of Car OB 1 . 

Second, the visual extinction toward the stars in Car OB 1 
is low, Av <3 (Humphreys 1978; Hoekzema et al. 1992). At 
this extinction and distance, a massive star should be clearly 
detected at the center of GSH 288.3-0.5-28. For exam- 
ple, the 06V star considered above would have a magnitude 
y < 9. We have performed an exhaustive search for such 
sources (e.g., [H umphreys 1978; Forte & Orsatti 1981; ReeJ 
'1998li'Kaltcheva 2003), from which the only catalogued mas- 
sive star projected near the center of GSH 288.3-0.5-28 is 
IX Car (HD 94096; see Fig. [TJ, an M2Iab stai- with a mass 
of ~ 20 M0. The distance estimated to I X Car of ~ 1 .6 kpc 
and its LSR radial velocity of -17 km s~' dHumphrevsl 1 978t) 
suggestthat it is a foreground object. Other massive stars in 
die vicinity, such as HD 93843, HD 94230, HD 305599 and 
LS 1976, afl lie on the very edge or outside of GSH 288.3- 
0.5-28 in projection. 

Finally, the photoionization from stars in Car OBI, most 
notably the 05III star HD 93843, should fully ionize the shell. 
If we assume that the shell thickness is 20% of its radius, 
the edges of the shell should show Ha emission at a sur- 
face brightness of thousands of rayleig hs. We have exam- 
ined sensitive image s of this region re.g.. lBuxton et al.lll998t 
■Gaustad et al. 2001); these show no Ha emission associated 
with GSH 288.3-0.5-28 down to much lower surface bright- 
ness limits. We therefore beheve that GSH 288.3-0.5-28 
is not exposed to the high ionizing fluxes associated with 
Car OBI, but likely has a higher extinction than and is be- 
hind these coincident stars. 

An important point is that the high extinction that we have 
invoked only need imply a slightly larger distance for the H I 
shell than for Car OBI, since extinction in this direction has 
little correspondence with distance. This is clear from con- 
sideration of the ^ 30 Wolf-Rayet (WR) stars within 3° of 
GSH 288.3-0.5-28 (van der Hucht 2001). For stars from this 
sample with distances in the range 1.7-3.7 kpc as appropriate 
for GSH 288.3-0.5-28, can be as high as 6-7. 

3.3. The Central Source Associated with GSH 288.3-0.5-28 

If GSH 288.3-0.5-28 is not associated with any of the 
known bright stars in this direction, what is powering it? If the 
extinction to GSH 288.3-0.5-28 is significantly higher than to 
the stars in Car OBI, its associated star might not be obvious. 
We first consider K or M supergiants, which should be promi- 
nent in the near-infrared — we estimate a 1.3-/im magnitude 
J < 6.5 for a red supergiant at the distance to GSH 288.3- 
0.5-28. Howeve r, using the 2MASS point source catalog 
(lCutrietal.ll2()0l . we find no candidates for such a source 
inside GSH 288.3-0.5-28, other than IX Car (see E3- We 
searche d for faint O B stars by considering the Tycho-2 cat- 
alog (H0g et al. 2000), which is 90% complete down to a V 
magnitude of 11.5. We assume that any massive star with 
4 in this well-studied region would have been classified 
in earlier efforts, and so only consider stars with Ay > 4. For 
standard reddening, this implies a color excess Eb-v ^1-3 
and hence for OB stars B-V > 1 . Just nine such stars from 
the Tycho-2 catalog lie within the perimeter of GSH 288.3- 
0.5-28, of which three are K or M stars with low extinction, 
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four have high proper motions which indicate that they are 
foreground objects, one is a foreground BIV star (LS 1956) 
which is heavily reddened, and one has near infrared mag- 
nitudes from 2MASS which are inconsistent with it being a 
massive star at the distance to GSH 288.3-0.5-28. 

It is certainly possible that stars too faint to be in the 
Tycho-2 catalog could also be candidates. While we thus 
accept that we cannot make a definitive identification using 
available data, there is one remaining source, namely the 
AXP IE 1048.1-5937, which we argue below represents a 
viable association with GSH 288.3-0.5-28. If IE 1048.1- 
5937 generates a wind powered by its spin-down, as do many 
radio pulsars, then its spin parameters imply only L„ ~ 4 x 
lO-'-' ergs s"' (Gavriil & Kaspi 2004). This is far too low to be 
responsible for the surrounding shell. However, an intriguing 
alternative is that GSH 288.3-0.5-28 was blown by the star 
whose collapse formed IE 1048.1-5937. As for associations 
between pulsars and SNRs, the likelihood of this association 
should be judged by geometric correspondence, agreement in 
distance and age, and evidence for physical interaction. 

First considering geometry, it is clear from Figure [2 that 
IE 1048.1-5937 is close to the center of GSH 288.3-0.5- 
28 in projection. Isolated neutron stars generally have high 
space velocities, > 300 km s"'( Arzoumanian et al. 2002),*' 
so that any neutron star older than ^100 kyr will have 
moved far from its progenitor's wind bubble. Thus the as- 
sociation is only viable if we can argue independently that 
IE 1048. 1-5937 is extremely young. Indeed there is good ev- 
idence that this is the case, since associations with SNRs for 
other AXPs argue that these sources all have ages < 10 kyr 
(^jaensler et al. 1999, 2001). In this case the AXP should 
be centrally located in any progenitor wind bubble, as ob- 
served. We note that the "characteristic age" of IE 1048.1- 
5937, f /2P « 4 kyr (Gavriil & Kaspil l2004). is indeed small 
(where P is the star's spin period).' However, we caution that 
this and other magnetars do not spin down like normal radio 
pulsars ( Thompson et al. 2000; Gavriil & Kaspi 2004), so this 
age estimate does not provide independent evidence for youth. 

Second, we consider distance estimates. The systemic ve- 
locity of GSH 288.3-0.5-28 puts it at 2.7 ± 1 kpc. No dis- 
tance estimates exist for IE 1048.1-5937, except for the fact 
that the absorbing column of hydrogen inferred from its X- ray 
specti-um, A^H « 1-1 X 10^^ cm"^ (e.g.. lTiengo et alJl2002h . is 
much higher than that for the Carina complex (which includes 
Car OBI) at a distance of k, 2.5 kpc. It has thus been argued 
that 2.5 kpc is a firm lower limit on the distance to this source 
iOzel et al. 2001). However, the high column to IE 1048.1- 
5937 is still compatible with the distance to GSH 288.3-0.5- 
28. For standard gas to dust ratios, A ^^ 1.1 x 10^^ cm~^ 
implies a visual extinction Ay « 5.8 (Wang & Chakrabartv 
12002). The WR catalogue of vanderHucht (2001) demon- 
strates that this level of extinction is consistent with distances 
in the range 2- 12 kpc in this region, so is easily reconciled 
with the distance to GSH 288.3-0.5-28. 

Third, we need to determine if the age of GSH 288.3- 
0.5-28 is consistent with that expected for the progenitor of 
IE 1048. 1-5937. A stai- with a ZAMS mass of ^30-40 Up, 
and so lar metallicity lives for 5-6 Myr (e.g. JMevnet et al.1 
(1991 . Since the age of the AXP is negligible in compari- 

* We note that no veloc ity measurements yet exist specifically for magne- 
tars IGaensler et all2001l> . 

' IE 1048.1-5937 experiences significant torque variations; we have 
adopted the long-term average for P. 



son, this sho uld a lso be the age of the surrounding bubble. 
However, in ^3.11 we esimated an age for the shell of 1 . 1 ± 
0.5 Myr There is clearly a large discrepancy between these 
two estimates. However, this same age problem has been ob- 
served for many other H 1 shells around individual massive 
stars (e.g., Gervais & St-Louis I_9S2i Cap pa_&_He rbstmeie3 
120001) . There is no simple resolution to this discrepancy, al- 
though in some cases it may be explained by a "blow out" 
into a lower density environment ( Oev & Smedlev 1998). We 
conclude that any discrepancy between the ages estimated for 
GSH 288.3-0.5-28 and for the progenitor of IE 1048.1-5937 
is not a strong argument against their association. 

Finally, we consider direct physical evidence that 
IE 1048.1-5937 and GSH 288.3-0.5-28 are associated. 
Young neutron stars should be embedded in young SNRs, but 
for IE 1048.1-5937 no such SNR is observed (Gaensler et aQ 
120011) . A simple explanation is that the associated SNR is 
expanding into a low density bubble, so that it does not yet 
produce observable emission (e.g., Ciotti & D'Ercole 1989). 
Thus from the absence of a SNR, it is reasonable to expect 
that IE 1048.1-5937 should show evidence for a surrounding 
cavity: GSH 288.3-0.5-28 clearly fulfills this prediction. For 
a SNR shock velocity of 5000 km s"\ the blast wave should 
impact the shell walls ^ 3 kyr after core collapse. The lack 
of any radio or X-ray emission from this event requires the 
neutron star to be younger than this, consistent with the small 
ages expected for AXPs. 

To summarize, a very young neutron star with a massive 
progenitor should be centrally located in an expanding Hi 
shell, with no evidence for a surrounding SNR. This is ex- 
actly what we observe for GSH 288.3-0.5-28 and IE 1048.1- 
5937. Both sources are consistent with being at a distance of 
~ 3 kpc, behind ^ 6 magnitudes of optical extinction. 

4. IMPLICATION FOR MAGNETARS 

As explained in ^ a fundamentally unresolved issue in 
the study of compact objects is why some neutron stars are 
ordinary radio pulsars, while others are magnetars. Since 
the initial mass function (IMF) sharply declines with increas- 
ing ZAMS mass, most neutron star progenitors will have 
masses near the minimum mass for core collapse, i.e., 8- 
9 Mq. We have presented evidence here that the progeni- 
tor of AXP IE 1048.1-5937 was considerably more massive 
than this. For some SGRs, possible associations with massive 
star clusters similarly argue for high-mass progenitors (e.g., 
Klose et al. 2004; Eikenberrv et al. 2004). 

We thus propose that the difference between normal pul- 
sars and magnetars is the progenitor mass. We note that mas- 
sive (> 25 M0) stars do not always form black holes: for 
solar metallicity, mass loss causes single stars heavier than 
~ 60 M0 to form neutron stars (Hegeret al. 2003); a pro- 
genitor of lower mass (25 - 40 Mq) whose binary companion 
strips its out er envelope before core c ollapse can also form a 
neutron star (iFrver & Kalogeral2001l) . 

Why should massive stars form magnetars? 
Duncan & Thompson ( 1992) and [Thompson & Dunca^ 
d 19931) argue that magnetars result from rapidly rotating 
(f ~ 1 ms) proto-neutron stars, in which an efficient large- 
scale dynamo operates in the first few seconds after birth, 
generating a super-strong magnetic field with significant 
multipolar components. Hegeret al. (2005) have recently 
carried out a series of calculations of differentially rotating 
magnetized supernova progenitors. They find that magnetic 
torques are especially effective at transferring angular mo- 
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mentum away from the stellar core during the red supergiant 
and helium burning phases of evolution. For progenitors 
with masses ^ 10-15 Mq, this results in neutron stars 
with initial periods ^ 10-15 ms, too slow to generate 
magnetar-like fields. However for more massive stars, the 
reduced time interval between hydrogen depletion and the 
supernova results in limited braking of the core, producing 
much more rapidly spinning neutron stars. For example, a 
35 Mq progenitor results in a neutron star of initial period 
3 ms, in the range needed to give birth to a magnetar. 

Within this scenario, we can immediately make a predic- 
tion as to the relative birth rates of radio pulsars versus mag- 
netars. For argument's sake, we suppose that the mass cut 
between normal pulsars and magnetars is at a ZAMS mass of 
25 Mq. Then for an IMF of slope a = 2.35 for massive stars 
(dN/dM cx M""; Massev 2003), only - 20% of progenitors 
can potentially form magnetars. Since ~ 50% of progenitor s 
heavier than 25 M© will form black holes jHeger et alJ2003l) . 
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Fig. 1 . — H I in the field of the AXP IE 1048.1-5937, from SGPS data. The three panels plot the surface brightness of H I as a function of Galactic longitude, 
Galactic latitude and LSR velocity. Each greyscale ranges between -10 and +78 Jy beam ' , as shown by the scale bar on the upper right. The spatial resolution 
of the data are 142" x 125", while the velocity resolution is 0.82 km s~^ The upper right panel shows the image plane at an LSR velocity of -28 km s"'; the 
positions of IE 1048.1-5937 and of the star DC Car are marked by a cross and a star, respectively. The left panel shows ab-v diagram at / = 288.°25; the lower 
panel shows al-v diagram at & = -0.°52. 



